Generalized Hartree-Fock Theory for Interacting Fermions in Lattices: Numerical Methods

We present numerical methods to solve the Generalized Hartree-Fock theory for fermionic systems in lattices, both in thermal equilibrium and out of equilibrium. Specifically, we show how to determine the covariance matrix corresponding to the Fermionic Gaussian state that optimally approximates the quantum state of the fermions. The methods apply to relatively large systems, since their complexity only scales quadratically with the number of lattice sites. Moreover, they are specially suited to describe inhomogenous systems, as those typically found in recent experiments with atoms in optical lattices, at least in the weak interaction regime. As a benchmark, we have applied them to the two-dimensional Hubbard model on a 10x10 lattice with and without an external confinement.Comment: 16 pages, 22 figure

Tensor network states and geometry

Tensor network states are used to approximate ground states of local Hamiltonians on a lattice in D spatial dimensions. Different types of tensor network states can be seen to generate different geometries. Matrix product states (MPS) in D=1 dimensions, as well as projected entangled pair states (PEPS) in D>1 dimensions, reproduce the D-dimensional physical geometry of the lattice model; in contrast, the multi-scale entanglement renormalization ansatz (MERA) generates a (D+1)-dimensional holographic geometry. Here we focus on homogeneous tensor networks, where all the tensors in the network are copies of the same tensor, and argue that certain structural properties of the resulting many-body states are preconditioned by the geometry of the tensor network and are therefore largely independent of the choice of variational parameters. Indeed, the asymptotic decay of correlations in homogeneous MPS and MERA for D=1 systems is seen to be determined by the structure of geodesics in the physical and holographic geometries, respectively; whereas the asymptotic scaling of entanglement entropy is seen to always obey a simple boundary law -- that is, again in the relevant geometry. This geometrical interpretation offers a simple and unifying framework to understand the structural properties of, and helps clarify the relation between, different tensor network states. In addition, it has recently motivated the branching MERA, a generalization of the MERA capable of reproducing violations of the entropic boundary law in D>1 dimensions.Comment: 18 pages, 18 figure

Entanglement renormalization and boundary critical phenomena

The multiscale entanglement renormalization ansatz is applied to the study of boundary critical phenomena. We compute averages of local operators as a function of the distance from the boundary and the surface contribution to the ground state energy. Furthermore, assuming a uniform tensor structure, we show that the multiscale entanglement renormalization ansatz implies an exact relation between bulk and boundary critical exponents known to exist for boundary critical systems.Comment: 6 pages, 4 figures; for a related work see arXiv:0912.164

Can One Trust Quantum Simulators?

Various fundamental phenomena of strongly-correlated quantum systems such as high-$T_c$ superconductivity, the fractional quantum-Hall effect, and quark confinement are still awaiting a universally accepted explanation. The main obstacle is the computational complexity of solving even the most simplified theoretical models that are designed to capture the relevant quantum correlations of the many-body system of interest. In his seminal 1982 paper [Int. J. Theor. Phys. 21, 467], Richard Feynman suggested that such models might be solved by "simulation" with a new type of computer whose constituent parts are effectively governed by a desired quantum many-body dynamics. Measurements on this engineered machine, now known as a "quantum simulator," would reveal some unknown or difficult to compute properties of a model of interest. We argue that a useful quantum simulator must satisfy four conditions: relevance, controllability, reliability, and efficiency. We review the current state of the art of digital and analog quantum simulators. Whereas so far the majority of the focus, both theoretically and experimentally, has been on controllability of relevant models, we emphasize here the need for a careful analysis of reliability and efficiency in the presence of imperfections. We discuss how disorder and noise can impact these conditions, and illustrate our concerns with novel numerical simulations of a paradigmatic example: a disordered quantum spin chain governed by the Ising model in a transverse magnetic field. We find that disorder can decrease the reliability of an analog quantum simulator of this model, although large errors in local observables are introduced only for strong levels of disorder. We conclude that the answer to the question "Can we trust quantum simulators?" is... to some extent.Comment: 20 pages. Minor changes with respect to version 2 (some additional explanations, added references...

Many body physics from a quantum information perspective

The quantum information approach to many body physics has been very successful in giving new insight and novel numerical methods. In these lecture notes we take a vertical view of the subject, starting from general concepts and at each step delving into applications or consequences of a particular topic. We first review some general quantum information concepts like entanglement and entanglement measures, which leads us to entanglement area laws. We then continue with one of the most famous examples of area-law abiding states: matrix product states, and tensor product states in general. Of these, we choose one example (classical superposition states) to introduce recent developments on a novel quantum many body approach: quantum kinetic Ising models. We conclude with a brief outlook of the field.Comment: Lectures from the Les Houches School on "Modern theories of correlated electron systems". Improved version new references adde